Magneto-optical and optical probes of gapped ground states of bilayer graphene 



(N 

o 

bJE> 
< 



G 

o 

> 

m 

(N 

O 

(N 



X 



E.V. Gorbar, 1,2 V.P. Gusynin, 2 A.B. Kuzmcnko, 3 and S.G. Sharapov 2 

department of Physics, Taras Shevchenko National Kiev University, 
6 Academician Glushkov ave., Kiev 03680, Ukraine 
A Bogolyubov Institute for Theoretical Physics, National Academy of Science of Ukraine, 
14~b Metrologicheskaya Street, Kiev, 03680, Ukraine 
3 DPMC, University of Geneva, 1211 Geneva 4, Switzerland 
(Dated: March 3, 2013) 

We study the influence of different kinds of gaps in a quasiparticle spectrum on longitudinal and 
transverse optical conductivities of bilayer graphene. An exact analytical expression for magneto- 
optical conductivity is derived using a low-energy two-band Hamiltonian. We consider how the layer 
asymmetry gap caused by a bias electric field and a time-reversal symmetry breaking gap affect the 
absorption lines. The limit of zero magnetic field is then analyzed for an arbitrary carrier density in 
the two-band model. For a neutral bilayer graphene, the optical Hall and longitudinal conductivities 
are calculated exactly in the four-band model with four different gaps and zero magnetic field. It 
is shown that two different time-reversal symmetry breaking states can be distinguished by analyz- 
ing the dependence of the optical Hall conductivity on the energy of photon. These time-reversal 
symmetry breaking states are expected to be observed experimentally via optical polarization ro- 
tation either in the Faraday or Kerr effects. We analyze a possibility of such an experiment for a 
free-standing graphene, graphene on a thick substrate, and graphene on a double-layer substrate. 

PACS numbers: 78.67.Wj, 78.20.Ls, 81.05.ue 



I. INTRODUCTION 

Since its experimental discovery^ bilayer graphene be- 
came a separate subject of research due to the fea- 
tures which make it unique among the other two- 
dimensional condensed matter systems. Its low-energy 
electron spectrum^ combines characteristics of monolayer 
graphene and traditional two-dimensional electron sys- 
tems. It consists of two inequivalent pairs of parabolic 
valence and conductance bands, touching each other at 
K and K' points, and charge carriers are massive and 
possess a chirality. 

A unique feature of bilayer graphene is that an electric 
field E±_ applied perpendicular to the layers results in the 
opening of a tunable gap between the valence and con- 
duction bands. The value of E± can be controlled exter- 
nally by chemical doping and gating. Since, in contrast to 
single layer graphene, the density of states remains finite 
even in the unbiased and undoped (neutral) bilayer, there 
are theoretical predictions^— that the electron-electron 
interaction can result in spontaneous symmetry break- 
ing and opening a gap even in the absence of a magnetic 
field. The nature of the gapped state is much debated 
in the literature. Possible scenarios include anomalous 
quantum Hall (QAH), quantum spin Hall (QSH), layer 
antifcrromagnet (LAF) states, etc. (see Ref. M for a gen- 
eral discussion). Technically, all these gapped states dif- 
fer in the way how they break an approximate 5*1/(4) 
spin-valley symmetry of the low energy Hamiltonian of 
bilayer graphene. 

In an external magnetic field symmetric bilayer 
graphene exhibits anomalous quantum Hall (QH) effect^ 
with the filling factors v — ±4n with n = 1,2, .... As 
in single-layer graphene, this QH effect is caused by the 



anomaly of the zero-energy lowest Landau level which 
is eightfold degenerate. The subsequent experiments^— 
showed that in higher magnetic fields the degeneracy of 
the lowest Landau level is completely resolved and new 
QH states with filling factors v — 0,±1,±2, ±3 appear. 
It turned out that the activation energy gaps for these 
QH states depend linearly on the magnetic field B. This 
behavior can be contrasted with the case of single layer 
graphene, where the corresponding gaps scale as \f~B~. 
As suggested in Refs.lMrTl. the difference between bilayer 
and single layer graphene is caused by a strong screening 
of the Coulomb interaction in the former. 

Interestingly, the experiment o 1 1 ' 1 5 demonstrated that 
neutral bilayer graphene remains gapped even when 
E± = and an applied magnetic field B vanishes. In 
particular, the authors of Ref. [TH concluded that the re- 
sults of their measurements are most consistent with the 
LAF state. Nevertheless, further theoretical and exper- 
imental work is necessary to ascertain the nature of the 
gapped state. 

Optical spectroscopy studies proved to be a useful tool 
in investigation of carbon systems (see Refs. [l8lfl9l for 
an overview). The rclativistic-likc gaplcss dispersion of 
quasiparticles in single layer graphene results in a uni- 
versal and constant optical conductivity <j xx = e 2 /4h 
for the photon energies above the threshold which is 
twice the Fermi energyi 2 ^— Magneto-optical properties 
of single layer graphene reveal themselves in the spec- 
troscopy of Landau levels transition a 29 ' 30 and the giant 
Faraday effect.— Optical properties of gapped single layer 
graphene were intensively studied also in a series of works 
both in zero and finite magnetic fieldi^— 

Optical methods turned out to be especially fruitful 
for bilayer graphene. Theoretical description of its op- 
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tical and magneto-optical properties 3 - 7 . - — involves more 
parameters as compared to single layer graphene and all 
these parameters were found experimentally^— in the 
B = case as well as in the presence of magnetic field 4£ 



As suggested in Ref. H^, optical methods could also 
be used for investigating the symmetry breaking gapped 
states in bilaycr graphene discussed above. A particu- 
lar case of the QAH stated which breaks time-reversal 
symmetry explores the idea that such state would show 
up in the rotation of the polarization of light. The pa- 
pers on magneto-optical conductivit y 37 ! 38 do not take 
into account the presence of gapped states in bilaycr 
graphene, while the paper— considers only the case of 
time-reversal symmetry breaking gapped state in zero 
magnetic field. In a recent paper ^ it was shown that 
the infrared and far-infrared absorption spectroscopy in 
bilayer graphene at zero or finite doping in zero magnetic 
field can distinguish gapped states from the gapless un- 
perturbed and ncmatic states due to their qualitatively 
different lineshapes. Since now there is an interest in 
crossover between various gapped states at finite and zero 
magnetic field, in the present paper we derive analytical 
expressions for magneto-optical conductivity which in- 
clude both arbitrary gapped states and finite magnetic 
field. 



The paper is organized as follows. In Sec. UH we in- 
troduce the model 4x4 and 2x2 Hamiltonians and 
discuss various types of gaps. In Sec. [ill] analytical ex- 
pressions for optical conductivities in a magnetic field for 
circularly polarized light are derived using the two-band 
model. We discuss how the opening of different kinds 
of gaps influences the conductivities. Then, the limit of 
zero magnetic field is analyzed for an arbitrary carrier 
density. In Sec. IIV1 we derive the optical Hall conductiv- 
ity in zero magnetic field in the presence of time-reversal 
symmetry breaking gap. In contrast to Ref. l48l . we use 
a four-band model to derive an exact expression for the 
transverse conductivity which allows us to probe a wider 
range of frequencies. We consider the general case where 
both time-reversal-invariant and -noninvariant gaps are 
present. The time-reversal-invariant gap can be gener- 
ated due to a perpendicular bias electric field, therefore, 
we study the optical Hall conductivity dependence on ap- 
plied bias electric field. In Sec. [V] we consider the detec- 
tion of the time-reversal symmetry breaking states in the 
optical polarization rotation experiment, viz., by measur- 
ing the Faraday or Kerr effects. We analyze a possibil- 
ity of such an experiment for a free-standing graphene, 
graphene on a thick substrate, and graphene on a dielec- 
tric layer on top of a thick layer. Finally, in the Dis- 
cussion section, we give a brief summary of our results, 
and the Appendix at the end of the paper contains an 
analysis of the optical spectral weight. 



II. MODELS AND NOTATION 

We consider bilayer graphene in the continuum approx- 
imation using both a 4 x 4 Hamiltonian in the absence of 
magnetic field B — and its low-energy 2x2 approxima- 
tion in the case B ^ 0. The effective Hamiltonian of the 
four-component model at fixed spin, s = ±, and valley, 
£ = ±, reads as 
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(1) 



where it = p x + ip yi vf ~ 10 6 m/s is the in-plane Fermi 
velocity, 71 = 0.38 eV is the inter-layer hopping between 
pairs of orbitals that lie directly below and above each 
other. In our study, we neglect the tight-binding param- 
eters 73 (which leads to a trigonal warping of the band 
structure at low energies) and 74 whose effects fall be- 
yond the scope of the present paper. 

The Hamiltonian acts on wave function with four com- 
ponents corresponding to the atomic sites Al, B2, A2, 
Bl in the valley K (£ = +1) and B2, Al, Bl, A2 in the 
valley K' (£ = — 1). In the effective Hamiltonian we in- 
cluded also quasiparticlc gaps A^ s dynamically generated 
due to the electron-electron interaction. 

The possible gapped states ^ 17 ' 50 can be classified con- 
sidering how the gap A^ s depends on valley and spin 
indices. Obviously, selecting gapped states symmetric 
and antisymmetric in valley and spin, the most general 
expression for the gap is given by 



A es = U + sU T + £A T + £sA, 



(2) 



where while gaps U and A are invariant with respect to 
the time reversal symmetry, gaps At and Ut are not. 

Among these four gaps the first gap U besides being 
dynamically generated can be induced by gating which 
creates a perpendicular electric field E±. In this case 
we will denote it as Uo(= e.E±d/2), where d = 0.35 nm is 
the distance between layers. Thus the gap U is related to 
the layer-polarized state caused by a potential difference 
between layers. Note that since the gap term in Eq. JT]) 
is introduced in the left upper corner, strictly speaking 
the gap U only approximates the layer asymmetry gap 
[see the discussion below Eq. (|3J)]. In the notations of 
Ref. H, the state related to the gap U is called quantum 
valley Hall (QVH) state. The gap Ut corresponds to 
the layer antiferromagnet state, At corresponds to the 
quantum anomalous Hall state and A is related to 
the spin-polarized or quantum spin Hall stated The cor- 
responding order parameters are summarized in Table [U 

All these order parameters except the first one were 
suggested for describing the ground state of bilayer 
graphene at the neutral point in the absence of external 
electric and magnetic fields. Indeed, while for a suffi- 
ciently large top-bottom voltage difference Uq, the layer 
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TABLE I: Possible gapped states in bilayer graphene at neutral point. QVH, quantum valley Hall; LAF, layer antiferromagnet; 
QAH, quantum anomalous Hall; QSH, quantum spin Hall. The summation over s — ± is implied. 



polarized QVH state is realized, the experiments^^ 
demonstrated that as the value Uo decreases there is a 
phase transition to another state. This eliminates the 
QVH state as a possible candidate for the ground state 
of bilayer graphene at the neutral point in the absence of 
external fields. 

On the other hand, for sufficiently large magnetic field 
B, the spin-polarized QSH state is realized. Recent ex- 
perimental data^ show the absence of a phase transition 
as B decreases to zero. This suggests that the QSH state 
could be the ground state of bilayer graphene at the neu- 
tral point in the absence of external fields. According 
to Ref. the LAF state is adiabatically connected to 
the QSH state at high magnetic field, therefore, the LAF 
state could be also the ground state of neutral bilayer 
graphene in the absence of external fields. 

To investigate the role of different kinds of gaps in the 
magneto-transport properties of bilayer graphene, it is 
convenient to use the effective low-energy Hamiltonian 
which was derived in Ref. fusing Green's functions. The 
2x2 Hamiltonian can also be obtained from the 4 x 
4 Hamiltonian by integrating out Bl, A2 fields which 
correspond to the Bcrnal stacked orbitals and is valid 
within the energy range |e| < 7i/4. The trigonal warping 
term neglected in Eq. (JTJ) also restricts the validity of 
the effective Hamiltonian at low energies^ In an external 
magnetic field this Hamiltonian takes the form 



£diag([7, — U). One can check that the energy spec- 
trum corresponding to the Hamiltonians Hu and H' v is 
practically identical. For zero magnetic field the com- 
bined effective Hamiltonian with the gap A^ s = U for 
two valleys, H(p, U) = H cB (£ = +l,p,U) ® E eS (£ = 
— l,p, U) is time-reversal invariant under the transfor- 
mation (IIi ® n)H*(p, U)(Ui <g) n) = H{-p, U), where 
ni swaps £ = +1 and £ = —1 in valley spaced As 
in the 4x4 case the presence of the gap A^ s = £Ay 
breaks the time-reversal symmetry. In the presence 
of a magnetic field the corresponding symmetry trans- 
formation becomes (Hi (2D Ti)H* (p, At, B)(Hi <E> t\) = 
H(-p,-A T ,-B). 

We calculate the optical conductivity analytically us- 
ing the Kubo formula, 



ftpg(n + io)-ng(o)] 



in 



(4) 



where n^(f2 + iO) is the retarded current-current cor- 
relation function obtained by analytical continuation 
[11^(0) = Hij(iil m — > fi + ie)] from its imaginary time 
expression, and fi is the energy of photon. Neglecting the 
vertex corrections, the calculation of the current-current 
correlation function reduces to the evaluation of the bub- 
ble diagram 



H 



eff 



2m 



(3) 



where tt — h(—iD x — D y ) is now expressed via the co- 
variant derivatives Di = di + (ie/hc)Ai with the electron 
charge — e < 0, and the effective mass of the carriers 
m = -fi/(2vp). The external magnetic field B = V x A 
is applied perpendicular to the plane along the positive z 
axis. The Hamiltonian © acts on a wave function with 
two components corresponding to the atomic sites Al, 
B2 in the valley K (£ = +1) and B2. Al in the valley 
K> (£ = -1). 

It has to be noted that different gap terms in the ini- 
tial 4x4 Hamiltonian may result in the same effective 
2x2 Hamiltonian. For example, two different Hamil- 
tonians H[j = £diag([/, — U, 0, 0) which corresponds to 
Eq. flU with Aj s = U and Hu = £diag(U,-U,-U,U) 
which takes into account the asymmetry between on-site 
energies in the two layers, result in the same H v = 



IIij(ifi m ) = -^J dre^ T J d 2 rd 2 r' 



(5) 



x tr 



j 4 (r)G(r,r',r)^(r')G(r',r,-r) 



where ji(r) = —cdH/dA l is the electric current den- 
sity operator, G(r, r',r) is the electron Green's function 
(GF), V is the volume (area) of the system, (3 = 1/T is 
the inverse temperature, fl m = 27rm//3, and tr not only 
takes care of the 4 x 4 or 2 x 2 matrices, but also in- 
cludes summation over the valley and spin indices. In 
the presence of a magnetic field the GF is not transla- 
tional invariant and a special care should be taken in 
treating the translation noninvariant phase of the GF as 
done in Sec. IIIIl When the magnetic field is absent, the 
GF's is translation invariant and one can directly go from 
Eq. ([5]) to the frequency-momentum representation of the 
polarization operator as done in Sec. IIV1 
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III. TWO-BAND MODEL: OPTICAL 
CONDUCTIVITY IN AN EXTERNAL 
MAGNETIC FIELD 

In this section, we will consider the spin singlet gap 
A ^ = U + £ At ■ The energies of the Landau levels in the 
two-band model are 



£„ e = -£A 6 n = 0,l, 
Eant; = aM n (_, M n? = ^A 2 + uj 2 n(n - 1), n > 2, 



(6) 



where a = ±, lo c = heB/mc = h 2 /ml 2 is the cyclotron 
energy, and I = sj Hc/eB is the magnetic length. Using 
the values 71 = 0.38eV and vf = 1-02 x 10 6 m/s from 
Ref. |H (see also Ref. [H) one can estimate the effective 
mass of carriers in bilayer as m — 0.032m e , where m e is 
the electron mass. Accordingly, the cyclotron energy is 
equal to u c ps 0.116 meV(m e /m)B[T] = 3.62meVB[T]. 



A. Green's function and calculation of IL,j(f2) 



Here z = r 2 /(2l 2 ), L"(z) are associated Laguerre poly- 
nomials (by the definition L" 2 (z) = L" 1 (z) = 0), and 
(i is the chemical potential. For brevity of notation 
the subscript £ in M n £ is omitted in what follows, i.e. 
M n = M n j. Since the 2x2 Hamiltonian is quadratic 
in 7r and 7T , the electric current operator contains the 
derivatives 



e 
rn 



7T+ 
7T 



3v(*) = 



m \ in 



-in' 




(9) 



The phase factors in Eq. ([5]) cancel and we obtain at finite 
temperature 



In the two-band model © in the Landau gauge A = 
(0,Bx), the GF G(r, r',u>) in the mixed coordinate- 
frequency representation has a form 



G(r, r', oS) = cxp ( —i 



. (x + x'){y-y') 
2l 2 



G(r-r», (7) 



where the translation invariant part of the GF is repre- 
sented as a sum over the Landau levels^ 



1 / Mr'T ± 



n— — oo 



(10) 



tr 



ji{r, r')G(r, r', iui n )jj(r', r)G(r', r, iw„ + zfi m ) 



G(r,w) 



1 



:,/2 



E 



1 



t _ 2 Qz) (^ + M~^A c )£»(z) 



h 2 {x+iy)- j2 



(8) 



where u„ = 7r(2n + l)//3 and 



ix(r,r') 



M r > r ) = — 
m 



-id x - d y 

idx + dy+i *-*'+$-y') 

— d x + id. 

_ x—x'+i(y-y') 
V 2p 



d x + id,, - 



■ x-x -i(y-y') 





x-x'-i(y-y') 





(11) 



Noting that j,(r, r') = jj(r — r') depend on the difference 
of the coordinates r — r', one can integrate over r + r' in 
Eq. (flT)|) canceling the volume factor V in the denomina- 
tor 



n— — oo 



(12) 



tr 



ji(r)G(r,iu n )jj(r)G(r,iuJn + i^l m ) 



Writing Eq. (|T2"j) we also used that G(r, iuj n ) is even under 
the transformation r — > — r, while ji(r) is odd. 

The calculation of the optical conductivity follows 
closely the corresponding calculation done for monolayer 
graphene , 32 ' 33 so we directly proceed to the final ex- 
pression. The only difference is that we introduce a fi- 
nite Landau level width T at the very end of the cal- 
culation in the final expressions for the conductivities. 



5 



The final L = expression for the polarization operator 
Hg(O) = n* x (Q.) ± iU^ y {Q) takes the form 



fc=0 



y - 



%(AM w )-n F (A'M H2 ) 



A.A'=± 



A'M fc+2 - AA/ fc+ i ± (fi + iO) (13) 



AA'A| 



ftAA'£A 5 
M fe+1 M fe+ 2 I " M k+1 M k+2 



± 



where np(bj) = l/[cxp((w — fi)/T) + 1] is the Fermi dis- 
tribution function. It is clear that the last ~ Q term in 
II± may only be present if the time-reversal symmetry 
breaking gap, At ^ 0. 



B. Magneto-optical conductivity 

It is convenient to consider the optical conductivities 

<r±(n) = a xx (n) ± ia xy (n) = h\n£(si) - ng(o)]/(iO) 

which correspond to the opposite circular polarizations 
of light. After introducing a finite width of Landau levels 
in Eq. (fl~3|) . the final result for the complex conductivities 
reads as 



tr±(fi) 



fc=0 



53 [n F (\M k+1 ) - n F (A'M fc+2 ) 



A,A'=± 



(14) 



1 



1 - 



AA'Af 



A'A/ fc+2 - AM fc+ i \ M k+1 M k+2 



AA'£A ? 



± 



Q t A'Mfc+2 ± AM fe+ i + 2iT 



The expressions for conductivities in the case of more 
general gaps A^ s are obtained from Eq. (|14[) by replacing 
A 5 with A^ s and inserting the overall factor 1/2 and 
summing over the spin variable. Equation (|14p is one 
of the main results of this paper which generalizes the 
results of Refs. I37ll38l to a finite A^ case. For A^ = 0, 
there is a practically overall agreement between Eq. (fT4")l 
and the corresponding expressions from Ref. I37ll38l except 
to the intensity of the to 1 transition. To verify this 
issue, in Appendix [A] we consider the behavior of the 
spectral weight. We show both that our result agrees 
with the behavior of the spectral weight for B = and 
that the weight is conserved when the chemical potential 
moves from the region Mq < /j, < M 2 to the region fi > 
M 2 . 

It has to be stressed that rather simple analytical re- 
sult (|T4l was possible to obtain because we used the ef- 
fective 2x2 Hamiltonian. The 4x4 consideration is 



much more involve d 40 ' 41 . Moreover, the neglected tight- 
binding terms, including the trigonal warping, can only 
be treated numerically^ which is beyond the scope of the 
present work. 

The scheme of Landau levels and allowed transitions 
for bilayer can be found, for example, in the review^ so 
that we can go directly to the discussion of the behavior 
of magneto-conductivity. In Figs. [T] and [5] we show the 
results based on the computation of Eq. (fl"4"|) for the real 
part of <J±(fl) in units of exo = e 2 /(Ah) as a function of 

measured in meV. The highest value of H = 80 mev is 
less than 71 /4. In both figures we set B = 3 T, T = 5 K 
and T = 2mcV. 

The reference cas o 37 ' 38 A^ = is shown in Fig. [1] 
To scrutinize the curves in this figure, we provide the 
energies of the Landau levels, viz. for the n = 2 
level |-E ±2 | = M 2 « 15.3 meV, for the n = 3 level 
|_E ±3 | = M 3 m 26.6 meV and for the n = 4 level 
|£' ±4 | = M 4 fa 37.6 meV. The long dashed (red) curve 
is for <7_|_(Q) and the dash-dotted (black) curve is for 
cr_(0). Both curves are for fj, = lOmeV which is less 
than E 2 . Accordingly, in the long dashed (red) curve we 
observe the lines with the energy M 2 — A/0.1 ~ 15.3 meV, 
M 2 + M 3 w 41.9 meV and M 3 + Af 4 m 64.2 meV. They 
correspond to the transitions from the Landau levels with 
the energy £0,1 to E 2 , from E- 2 to E 3 , and from E—3 
to -E4, respectively. Since for the opposite polarization of 
light presented by the dash-dotted (black) curve there is 
no transition between the 0th (n — 0, 1) and E 2 level, the 
first observed line is for the transition from E^% to E 2 
levels with the energy M 2 + A/3 coinciding with the sec- 
ond line on the previous curve. The second observed line 
corresponds to the transition from i?_4 to E3 levels with 
the energy M3 + M4. The solid (blue) curve is for cr+(f2) 
and the short dashed (green) curve is for cr_(f2), and 
both curves are for \x = 20meV, so that E 2 < /x < E3. 
Accordingly, in the solid (blue) curve the transition with 
the energy M 2 — M\ is impossible, but instead the transi- 
tion with the energy M3 — M 2 w 11.3meV develops. Two 
remaining transitions from E^ 2 to E3, and from i?_3 to 
E4 which coincide with the transitions on the considered 
above long dashed (red) and dash-dotted (black) curves 
are also present. For the opposite polarization of light 
presented by the short dashed (green) curve the only pos- 
sible transition is from E^4 to E 3l so that the peak at 
the energy M3 + M4 is the only peak shared by all curves. 



C. Comparison of the QVH and QAH states 

Now we are at the position to discuss how the open- 
ing of either time-reversal symmetry breaking gap At or 
preserving this symmetry gap U changes the presented 
above picture (see Fig. (2). 

Again long dashed (red) curve is for o + (f2) and the 
dash-dotted (black) curve is for cr_(f2). Both curves are 
for = lOmeV, U = 5mcV and A T = 0. The solid 
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FIG. 1: (Color online) The real part of the conductivity a±(fl) 
in units of ao = e 2 /(4ft) versus the photon energy Q. in meV 
for magnetic field B — 3T, temperature T = 5K, and scat- 
tering rate T = 2meV. Long dashed (a + ) and dash-dotted 
(cr_) are for the chemical potential fi — lOmeV. The solid 
(o + ) and short dashed (<r_) are for the chemical potential 
H = 20meV. 



(blue) curve is for <7 + (f2) and the short dashed (green) 
curve is for <r_(f2), and while the value of the chemical 
potential is kept the same, fi = lOmeV, the gaps are now 
U = and At = 5meV. For relatively small values of 
the gaps, 5 meV, the energies of the levels with n ^ 0,1 
remain practically intact, viz. \E±2\ = Mi ss 16.1 meV, 
\E±z\ = M 3 « 27.1 meV and \E± 4 \ = M 4 ~ 37.9 meV. 
Accordingly, the last two peaks from the right remain 
practically unchanged, the first from the right peak with 
the energy M3 + M4 shifts in energy from 64.2 meV to 
65 meV, and the second from the right peak with the en- 
ergy M2+M3 shifts in energy from 41.9 meV to 43.2 meV. 
The most essential changes occur due to the partial re- 
moval of the degeneracy of the lowest Landau level. 

Since we took ji = lOmeV all curves in Fig. [2] have 
to be compared with the two curves in Fig. [T] plotted for 
the same value of the chemical potential, viz. the long 
dashed cr+(f2) curve and the dash-dotted <7_(fi) curve. 

Let us begin with the cr + (f2) polarization. We observe 
that for the gap U — 5meV the peak at 15.3 me V cor- 
responding to the transition from £0,1 to E2 levels is 
split into two peaks at 11.1 meV and 21.1 meV. These 
new peaks correspond to the transitions from the lev- 
els £ ,i = U and £"0,-1 = —U to £ 2 level, respec- 
tively. For the time-reversal symmetry breaking gap, 
At = 5 meV, we observe that there is only one peak at 
21.1 meV which corresponds to the transition from the 
only level £0.-1 = —At to £2. One can verify that the 
position of the peak is sensitive to the sign of the gap, 
i.e. for Ay = — 5meV the position of the peak is at 
11.1 meV, because there is the only level at the positive 
energy, £ ,i = -A T > 0. 

For the a-(Q) polarization, the discussed above transi- 
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FIG. 2: (Color online) The real part of the conductivity a±(Q) 
in units of ao — e 2 /(4ft) versus the photon energy Q in meV 
for magnetic field B — 3T, temperature T = 5K, scattering 
rate F = 2meV, and the chemical potential fj, = lOmeV. 
Long dashed (o + ) and dash-dotted (er_) are for the layer 
asymmetry gap U — 5meV and the QAH gap At = 0. The 
solid (er+) and short dashed (ff— ) are for the layer asymmetry 
gap U — and the QAH gap At = 5meV. 



tions are forbidden. The only small difference is observed 
between the height of the second from the right peaks in 
the solid (blue) and short dashed (green) curves. These 
curves interchange when the sign of the gap At is re- 
versed to — At- 

It should be noted that although a measurement of 
the optical Hall conductivity a xy (Cl) allows one to distin- 
guish holes and electrons,— some essential features such 
as splitting of the absorption peak into two peaks de- 
pending on the type of the gap remain present even in 
the diagonal conductivity a xx (Q). 

The obtained results allow us to conclude that the in- 
vestigation of the optical magneto-conductivity may pro- 
vide an additional insight into the nature of the gaps in 
the bilayer graphene. The features discussed here could 
be observed if the condition T < |A^| is satisfied. Since 
our value T = 2 meV is almost 5 times smaller than that 
used in Ref. |45| for fitting the data on the gated bilayer 
at B = 0, the observation of the discussed here effects 
requires samples of a much better quality. 



D. Zero-field and dc limits 

Starting from Eq. (fLlj) one can reproduce correctly 
B — > limit. Similarly to the consideration done for 
monolayer graphene^ introducing a continuum vari- 
able uj instead of M n given by Eq. and replacing 
the sum over n by the integral, in the limit T — > 
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for the real part of the diagonal optical conductivity 
&xx(ty = + (7 -(^))/2 one obtains 



Rea xx (n) 



ne 
~2h 



4Tcosh 2 ^ 



4Tcosh 2 ^ 



+ + ^k] sinh(|^|/2r)g(|0|-2|A^|) | 
2^ 9? J cosh(^/T) + cosh(0/2T) J 

(15) 

Here we also included the dependence of the gap Aj s 
on the spin variable s. The first term in curved brack- 
ets corresponds to the Drude peak while the second one 
describes interband electron-photon scattering processes. 
Finally, at T = we obtain 



Rea xx (tt) 



ire 



E 



ft 2 + 4A 2 
+ 40 2 es g(|0|-2max(| M [,[A^[) 

(16) 

This expression for longitudinal optical conductivity is 
similar to that obtained for single layer graphene 2 ^ and 
for topological insulators^ It is worth noting that there 
is a more deep analogy between the band structures and 
optical conductivities of the single layer graphene with 
Rashba term and biased bilayer graphene^ 

Similarly, one can consider the B — > limit for the 
optical Hall conductivity, a xy (fl) = (a + (Q)—<j-(Sl.))/(2i) 
which takes the form 



a xy (Q) 



4^ 
~h 



E 

Z,8=± 



I A* 



4lo 2 -{n + iO) 2 ' 



(17) 

Note that the last expression is an even function of chem- 
ical potential [i in contrast to the Hall conductivity in 
a magnetic field which is an odd function of /i. It is 
clear that a xy (fl) ^ only if the time- reversal symmetry 
breaking gaps Ut or Ay are nonzero. For zero tempera- 
ture we obtain 



>XIJ 



h(fl + iO) 



E t** 



x In 



£, a =± 

2max(|/z|, \A is \ 



+ fl + iO 



(18) 



2max(|/x|, \A 6s \) - (fi + zO)' 



Equation (fl8|) resembles the result obtained for topo- 
logical insulators.— Making the change of the variable, 
s — s- £s, in Eqs. ([TB|) and (|18l) . one can see that both con- 
ductivities are invariant under the interchange Ut A. 



The real part of the dc Hall conductivity takes the form 

A, 



Rea xy = -r E ^ 



h ^ ± ^max(|/x|, |A^|) 



h 



E^ 



sgn(A^), \n\ < |A Ss |, 
Ia*1 > |A 5s |, 



(19) 



which at fj, = 0, [/ = Ut = A = 0, and At ^ is in 
agreement with the corresponding expression in Ref. l58l . 
At the neutral point, /i = 0, we have, 



Rc a x 



ve 
~~h 



£sgn(A ?s ). (20) 



Clearly, the factor can take the values i/ = 0, ±2, ±4 de- 
pending on the relations between the gaps U, Ut, At, A 
and their signs. The QAH gap At and LAF gap 
Ut break time-reversal symmetry and this is the nec- 
essary condition for observation of nonzero Hall con- 
ductivity. For example, v = 4 is realized in the case 
A T > 0,U = U T = A = (for A T < 0, obviously, 
v = —4), thus QAH phase has a zero- field quantized 
charge Hall conductivity. In the QSH phase with A ^ 
and U = Ut = At = 0, two spin components have 
opposite Hall conductivity, hence zero charge Hall con- 
ductivity. On the other hand, in this phase a spin Hall 
conductivity <J xy {s = +) — u xy {s = — ) is nonzero and 
quantized. For At = and taking without the loss of 
generality all other gaps positive, we obtain v = —2 if 
the following conditions are satisfied 



U+ U T > A > \U-U T \, 
U + A > U T > \U- A\ 



(21) 



[v = 2 is obtained if we invert the signs of all three gaps 
and replace the gap values by their absolute values] . The 
second inequality in (|2ip arises as a consequence of the 
symmetry a xy (fl) under the interchange Ut 4-> A. 

In Fig. |3] we plotted the real part of frequency depen- 
dent Hall conductivity at zero temperature and scatter- 
ing rate 2r = 5 meV. The scattering rate T is introduced 
by replacing iO in Eq.(fl8]) by 2zT. In Fig. [3] two differ- 
ent cases for time-reversal breaking gaps are presented: 
1) Long dashed (red) line is for a finite QAH gap At = 
5mev, all the other gaps are zero, U = Ut = A = 0, 
and 2) For solid (blue) line the gaps are taken At = 0, 
U = 5mev, Ut = 4mev, A = 1.5 mev. Since the two- 
band model is valid for SI < 7i/4 ~ 0.1 eV, only this 
range of energies is shown in the figure. 

These two types of time-reversal symmetry breaking 
gaps in bilayer graphene are expected to manifest them- 
selves through a nonzero dc Hall response in the absence 
of a magnetic field. At finite frequencies, the behavior of 
the Hall conductivities differs essentially for two types of 
gaps: the curve for nonzero gap A crosses zero, while the 
curve for At =^ does not. The crossing takes place in 
infrared range of frequencies, while at near-infrared and 
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FIG. 3: (Color online) The real part of the optical Hall 
conductivity a xy (Q,) in units of e 2 /h versus the photon energy 
f2 in eV for T = fj, = 0, the scattering rate 2F = 5 meV. The 
long dashed (red) line is for the QAH gap At = 5mev, the 
other gaps are zero, U = Ut — A = 0. The solid line is for 
At = 0, the layer asymmetry gap U = 5mev, the LAF gap 
Ut = 4mev, and the QSH gap A = 1.5 mev. 



optical frequencies, ft ~ 7i/4 |Af s |, the Hall conduc- 
tivity is small and of order TLeo- xy ~ s A^ s \A^ s \/il 2 
[see, Eq. ([18])]. Nevertheless, it is possible to observe such 
gaps in a Hall response at optical frequencies as was sug- 
gested in Ref. HH. The point is that the Hall conductivity 
at optical frequencies is dominated by transitions to the 
high-energy bands and can be manifested in the polar 
Kerr or Faraday rotation. However, to make reliable cal- 
culations, it is necessary to move up to the four-band 
model. We study the optical Hall conductivity in the 
four-band model in Sec. IIVI 



IV. OPTICAL HALL CONDUCTIVITY IN THE 
FOUR-COMPONENT MODEL FOR ZERO 
MAGNETIC FIELD 

In this section, we will calculate and analyze the opti- 
cal Hall conductivity for bilayer graphene at the neutral 
point, i.e. [i = is set from the beginning. It was shown 
Ref. |48| that if At 7^ 0, the ac Hall conductivity, cr xy (H) 
exhibits a resonant enhancement at f2 = 71 which corre- 
sponds to the optical frequencies due to transitions from 
the low-energy bands to the high-energy bands. These 
optical interband transitions were effectively introduced 
into the two-band models using projector operators. 

Here, instead, we use the full four-band model ([1} and 
obtain an exact expression for optical Hall conductivity 
valid in a wide range of the photon energies. This enables 
us to investigate how behavior of a xy is affected by the 
other types of the gaps both in the infrared and optical 
ranges of frequencies. 

In the four-band model, the GF for quasiparticles at 
fixed spin and valley equals 



G(u),p) = (lj - H)- 



{u 2 - E 2 {p)){u 2 - E 2 {p)) 




(22) 



where 



2 2 



± 



(7i 2 -A? s ) 2 



(23) 



<p2( 7 -f + A|J, 



and the block matrices are 



J 



A 1 = 



(w + ^ s )( W 2 - 7 ?)- w 4p5 

VpllijPx +iPy) 2 



A 3 = £ 



v F li{Px -iPyf 



(w-£A es )(G 



-ll) 



■ ujv F p 2 



v f [uj(oj + £A Cs ) - v 2 F p 2 ]{p x - 
vf7i( w ~ CA ?s )(p x + ip y ] 



iPy) 



I uj?7i(w + £A £a )(p x - ip y ) 
y v F [w(u - CA^ S ) - v 2 F p 2 ](p x + ip y ) 

I vf7i( w + £,^s){Px + ipy) v f [uj(uj - £A ?S ) - v 2 F p 2 ]{p x - ipy 
y v f [lo(uj + £A ?s ) - v 2 F p 2 ](p x +ip y ) w_f7i(w - ^A ?s )(p x -ip y ), 

(u + ^){u(u - £A ?S ) - v 2 F p 2 ) 7l (w 2 
71 {to 2 



A? ^ 



A? 1 

is, 



{uj - £A ?s )( w (aj + ZA Ss ) - v 2 F p 2 ] 



(24) 



For |A^ S | -C v F p <C 71, we find for positive energies, 
£ l( p)^7 1 + ?M, ^(p)^ /(^p + ^Pl, 



i.e., Ei(p) and ^(p) describe the positive energy 
branches of the high and low energy bands, respectively 
(the negative energy branches of these bands are obtained 
multiplying the above expressions by —1). 
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According to the Kubo formula, the tensor of conduc- 
tivities is given by Eq. @, where at zero temperature 

/d d 2 
— ftr[j 1 G(w,p)j J G(w-n j p)] , 

(26) 

and the trace includes the summation over spin and 
valley degrees of freedom. Here, the current density 
j = edH/dp, or in components, 



jx = e&F 



/0 l\ / -i \ 

10 . . i 

10 9 0-i00 

\ 1 / \i / 

(27) 

Using Eq. (f2l)]) and taking the trace, we get 



n£,(fi) 



2e 2 v 2 F Q 



dujcPp N((jj, p) 



€=±,«=± 



)3 £>((„, p)' 

(28) 



where the numerator and denominator are, respectively, 



N(u, p) = v 4 F p A + v 2 F p 2 (4 7l 2 - lu 2 - {u - nf) 
+ uj 2 (uj - Q) 2 + A 2 s w(cj - O) 
+ 7i 2 (A| s -3w(o;-0)-0 2 ), 

D(w,p) = [w 2 - £ 2 (p)K - S|(p)] 

x [( W - O) 2 - i? 2 (p)][( W - fi) 2 - £ 2 2 (p)] . 



(29) 



Eq. ([28)1 implies that the function H^ y (O) is an odd func- 
tion of the energy: U^ y (-fl) = -U^ y (fl). On the other 
hand, from Eq. (H| we obtain that 11^.(0) = n^(-O), 
therefore, for the Hall conductivity we have the relation- 
ship a yx (Vt) = -a xy (£l). 



The integration over w and the angle in Eq. (|2"5)l can 
be done straightforwardly and we obtain the final general 
expression for the optical Hall conductivity in the four- 
component model of bilayer graphene, 



(T X y(Q) 



2h z -' 

f=±.s=± 



i 



El) 2 



1 



— ( 1 1 

£ 2 ^2^i - fl - 2iT + 2E 1 + n + 2iT 



E 2 \2E 2 -fl- 2iT 2E 2 + tt + 2iT 
4 7l 4 - 7l 2 (3£ 2 + 2E 1 E 2 + 3E 2 - 4A 2 J + (g - E 2 ) 2 - (E l - E 2 ) 2 A 2 S 



1 



2E 1 E 2 (E 2 - E 2 ) 2 
1 

Ei + E 2 - fl - 2iT + E 1 +E 2 + fl + 2iT 

I 



(30) 



Here as before we introduced a phenomenological impu- 
rity scattering rate T. In this form, the physical meaning 
of three terms in Eqs. (|30[) is quite transparent. Clearly, 
the first term in the square brackets describes transitions 
between the negative and positive branches of the high 
energy band Ei (p) , the second term is related to transi- 
tions between the negative and positive branches of the 
low energy band E 2 (p) 1 and the last term in curly brack- 
ets describes intcrband transitions. 



the form: 



where 



4e 2 z 



dy 



1 + M X M 2 



MiM 2 (Mi + M 2 f 



y- 



3MtM 2 + Ml + Ml + MfM. 



2 1 



M?M%(M 1 + M 2 ) 3 



(31) 



M lt2 (y,z)-- 



y± 



(i-z 2 f 



+ y(l + z 2 ).(32) 



The dc Hall conductivity in clean sample is obtained 
from Eq. (|30[) setting T = and £7 = 0. It can be rewrit- 
ten in terms of the dimcnsionless variables y = x/jf, 



At/7i (2?/ = Ml 



M 2 



1 — z ) and for the time- 



reversal symmetry breaking gap, A^ s = £Ar, it acquires 



To study the case At <C 71 we change y — 
then take the limit z — > 0. We find that a xy 



zy and 
= Ae 2 /h. 

The numerical study of Eq. (f3Tj) shows that a xy does not 
depend of z and always equals 4e 2 //i. Actually, this is 
a reflection of the fact that the dc Hall conductivity u xy 
can be written in terms of the topological Pontryagin 
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index 



C 



24tt 2 



dwd 2 ptv [Gd fl G~ 1 Gd u G~ 1 Gd\G~ 1 ] , 

(33) 



where d^G 1 = dG 1 /dp t _ l , p^ = (u,p x ,p y ), e^vX is the 
antisymmetric tensor, and G = (iui — H)~ l is the Green's diagonal optical conductivity 



function on the imaginary frequency axis^ 

In the case of a more general gap A^ s = U + £At the 
behavior of the dc Hall conductivity is different : it is zero 
for A T < U and a xy = 4e 2 /h for A T > U. 

We also provide the analytic expression for the complex 



E 



d.v 



ihSl^J (El -Elf 
t-± o 



2 7 fx 



El + Af 



1 



El \2E 1 -Q- 2iT 2E 1 + Q + 2iT 2E Y - 2iT 2E X + 2iT 



£| + A| 



1 



1 



1 



1 



2E 2 -£l- 2iY 2E 2 + n + 2iT 2E 2 - 2iT 2E 2 + 2iT 
1 1 1 



2EiE 2 \E 1 +E 2 -n-2iT E 1 +E 2 +Q, + 2iY E l + E 2 - 2iT Ei+E 2 + 2iT 
x [2E x E 2 (El - Etf-^lE 1 E 2 {E 2 1 + E 2 2 - 7l 2 ) 

+A 2 [(El + El)(E x - E 2 f - 5 7l 2 (£ 1 2 + E 2 2 ) + Q 1 2 1 E X E 2 + 4 7 2 ( 7 2 + A 2 )] -Af (Er - E 2 f] 



(34) 



Equation (j33J) is derived from the four-band model simi- 
larly to the optical Hall conductivity. We will use this in 
Sec. [V] for the analysis of Kerr and Faraday rotations. 



A. Optical Hall conductivity for the QAH state 

We immediately conclude from Eq. (|30p that for the 
QVH and LAF states, a xy (fl) = 0. For the QSH state, 
the optical Hall conductivity is also zero due to the sum- 
mation over spin in Eq. (|30[) . However, it is not equal 
to zero for the QAH state (A $s = £A T ). In Fig. H the 
dependence of Kca xy (U) is plotted. The thick curves 
are computed using Eq. (|30[) obtained in the four-band 
model. For comparison with Ref. we took the same 
values for the parameters 71 = 0.4 eV, 2T = 0.05 eV and 
plotted thin curves using the approximate expression (12) 
derived in Ref. |H. We also considered the two values of 
the gap: At = 1 meV (long-dashed (red) curve) and 
Ay = 0.5 meV (solid (blue) curve). We observe that the 
approximate expression is in agreement with the four- 
band model in the vicinity of Q = 71 and for higher en- 
ergies. For low-energies the behavior of the Hall conduc- 
tivity is correctly described only by the four-band model. 

In addition to the QAH gap At, the gating induces a 
finite time-reversal invariant gap Uq. Our multigap ex- 
pression (|3"0")) for the optical Hall conductivity also allows 
to study this case if one takes a gap A^ s = U + ^At- The 
results of the computation are presented in Fig. [SJ The 
QAH gap At is set to be 1 meV for all three curves, while 
the value of the layer asymmetry gap Uq is changed. The 
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FIG. 4: (Color online) The real part of the optical Hall con- 
ductivity a xy (Q.) in units of e 2 /h as a function of photon en- 
ergy fi. The thick lines are c omp uted using Eq. (|30|) and thin 
lines using Eq. (12) of Ref. |48|. The long-dashed and solid 
curves are for the QAH gap At = 1 meV and At = 0.5 meV, 
respectively. 



long dashed (red) curve is for U — Uq = 0.5 meV, the 
solid (blue) curve is for U = Uo = 5 meV, and the dash- 
dotted (black) curve is for U = Uq = 50meV. It is seen 
that in general the optical Hall conductivity is sensitive 
to the external electric field E± which decreases the value 
of the jump of a xy (£l) near f2 = 71 and shifts its position 
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FIG. 5: (Color online) The real part of the optical Hall con- 
ductivity a xy (Q) in units of e 2 /h as a function of photon en- 
ergy Q,. The QAH gap At = 1 meV. The long dashed curve 
is for the layer asymmetry gap U = 0.5 meV, the solid curve is 
for U = 5meV, and the dash-dotted curve is for U = 50meV. 



to the higher energies. 



The optical Hall conductivity for At 
superposition of three gaps 



and 



The Hall conductivity Rea xy (£l) is nonzero if the time- 
reversal symmetry is broken. This however is a necessary, 
but not sufficient condition. Indeed, although the LAF 
state gap Ut breaks time-reversal symmetry, it follows 
from Eq. (pRU that Re a xy (fl) — in this state. Then it 
seems that the Hall conductivity could be nonzero only 
when the QAH gap At 7^ 0. However, the analysis per- 
formed in Sec. IIIIDI shows that the dc Hall conductivity 
might still be nonzero even if A^ = 0. According to 
Eq. (|2ip . this happens if three other gaps Ut, U, and A 
are not zero and satisfy certain inequalities. 

As we mentioned in Sec.HJ the two states, LAF (A^ s = 
sUt) and QSH (A^ = £sA), are the most likely candi- 
dates for the gapped ground state of bilayer graphene at 
the neutral point in the absence of external fields. The 
gating of bilayer induces the layer asymmetry gap Uq. 
In the case Uo ^ 0, the analysis of the mean-held gap 
equations shows that the QSH state transforms into the 
superposition of the QSH gap A and the layer asymmetry 
gap U, while the LAF state into the the superposition of 
the LAF gap Ut and U, respectively. Furthermore, in the 
QSH state with two gaps a third LAF Ut gap may also 
open in the presence of magnetic impurities. Similarly, 
there may be a mechanism to generate a QSH gap, A for 
the LAF state with two gaps. Since these additional third 
gaps naturally should be much smaller than the other 
two, in this subsection we will consider the optical Hall 
conductivity for the two cases of superpositions of three 
gaps. In the first case, Ut is much smaller than A, U 



and, in the second, A is much smaller than Ut, U. How- 
ever, because of the symmetry of the Hall conductivity 
under the interchange Ut O A, the second case is in fact 
equivalent to the first one. In Fig.|5]we plotted the depen- 
dence Rea xy (SY) for fixed gaps Ut = 1 meV, A = 5 meV 
and for three different values of the layer asymmetry gap 
U = 5meV, lOmeV, 20meV and T = 25meV. We ob- 
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FIG. 6: (Color online) The real part of the optical Hall con- 
ductivity a xy (Q.) in units of e 2 /h as a function of photon en- 
ergy fi. The QAH gap A T = 0. The LAF gap U T = 1 meV, 
and the QSH gap A = 5meV. The long dashed (red) curve 
is for the layer asymmetry gap U = 5meV, the solid (blue) 
curve is for U = lOmeV, and the dash-dotted (black) curve 
is for U = 20meV. 

serve that the behavior of Kea xy (fl) becomes more rich 
than in pure QAH state and in addition to the jump at 
fl = 71 a new kink at lower energy develops. For smaller 
values of U = 5 — 10 meV this kink results in the two 
additional zeros on the function a xy (fl). 

This shows that if the Faraday or Kerr rotation is ob- 
served in zero magnetic held, not only can one state that 
some time-reversal symmetry state is present, but also 
shed a light on the specific nature of this state, by study- 
ing the dependence a xy (Q). To plot the figures in this 
section, we used a high value of the scattering rate T 
which shows that the observed features are rather robust 
with respect to the disorder. This allows us to conclude 
that these effects can be observed experimentally even if 
the values of the gaps are smaller than T. 



V. RELATION BETWEEN THE OPTICAL 
HALL CONDUCTIVITY AND THE 
FARAD AY/KERR ROTATION ANGLES 

Spontaneously broken time-reversal symmetry in bi- 
layer graphene should be observable via optical polariza- 
tion rotation when light is transmitted through the sam- 
ple (Faraday effect) or reflected by it (Kerr effect).— If 
graphene is deposited on a substrate, the rotation angles 
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depend not only on the optical conductivity of graphcnc, 
but also on the substrate properties. In this section, we 
provide formulas and calculate Faraday and Kerr rotation 
spectra for three practically relevant situations shown in 
the inset of Fig. [3 namely, (i) for free-standing graphene, 
(ii) for graphene on a thick substrate with a refractive 
index n = 1.5, (which closely matches the properties of 
SiC>2 and boron nitride in the spectral range of interest), 
and (hi) for graphene on a dielectric layer with n = 1.5 
and a thickness of d = 300 nm on top of a thick layer 
with n s = 3.5 (which corresponds to the most commonly 
used SiC>2/Si substrates). 

If graphene is deposited on a thick substrate with re- 
fractive index n, the experimental Kerr and Faraday ro- 
tation angles are defined, respectively, by the relations: 



-IK 



arg(r_) - arg(r + ) 



arg(t-) -arg(i+) 



where 



1 — n — Zq<j± 

+ Z <7± ' 



1 



t- 



1 



Z (T± 



(35) 



(36) 



(37) 



are the reflection and transmission coefficients at the 
'vacuum-film-substrate' interface, a± is the optical con- 
ductivity for the right ('+') and left ('-') circularly po- 
larized light calculated using Eqs. (|34|) and (f30|). and 
Zq = 4n/c is the impedance of vacuum. The case of 
free standing graphene is covered by the same relations 
if n is set to 1. 

In the limit Zo|<7±(f2)| <C n — 1 the formulas (|3"5|) and 
([3"6l are greatly simplified (see also Refl4§h: 



2Z Rca xy (n) 
n 2 -l ' 



(38) 



standing graphcnc, since the reflection coefficient itself is 
proportional to a 2 and therefore very small^ 

Next, we consider graphene on a double-layer sub- 
strate. In this case, Eqs. ([33]) and ([31)]) are still valid, 
but the reflection and transmission coefficients are calcu- 
lated differently: 



t± 



n>i± + 



toi±ri 2 T 2 t 10 ± 



1 - T 2 r 12 r 10 ± ' 
toi± T ti2 



(40) 



1 - T 2 r 12 r 10 ± 



where 



roi± = 
rio± = 
ri2 

T 



1 — n — Zqg± 
Z cr± ' 

Zq(7± 



I + 11 

n — 1 



toi± = 



1 + n - 
n — n s 



Z a± ' 

*12 = 



t 



1 + n + Z a± ' 
2n 



io± 
2/) 



1 + n + Z a± ' (41) 



n + n s n + n s 

cxp I c)nd} . 



Here, r,j-± and rji± are the reflection and transmission 
coefficients at the interface between media i and j (0, vac- 
uum; 1 and 2, the first and the second substrate layers) 
and r is the transmission coefficient for the first substrate 
layer. 

This calculation takes fully into account the Fabry- 
Perot interference in the 300 nm layer but not in the 
thick substrate. The corresponding results are shown by 
the short-dashed curves in Fig. [71 Due to the Fabry- Perot 
effect, the Faraday angle and especially the Kerr angle are 
no longer determined by Hea xy (fl) only. Interestingly, 
the Kerr angle above 0.4 eV is even inverted with respect 
to the case of a simple substrate. 



VI. DISCUSSION 



n+1 

from where it is obvious that both angles are proportional 
to the real part of a xy (Q,). Note that the approximation 
([35]) is not correct for the free-standing graphene. 

In Fig. [TJwe show the calculated Faraday and Kerr ro- 
tation spectra for free-standing graphene (the solid line) 
and for graphene on top of a thick substrate with n = 1.5 
(the dashed line), expressed in units of the fine structure 
constant a. The exact relations ([33]) and (|3S| were used. 
As an example, we take the case, where At = 1 meV and 
other gaps are equal to zero (Fig. 2]). One can see that the 
Faraday angle docs not differ much in the two cases and it 
indeed matches the real part of the Hall conductivity. On 
the contrary, the Kerr rotation for the free-standing sam- 
ple is about 10 times larger than for supported graphcnc. 
This does not necessarily mean, however, that the Kerr 
geometry is favorable to detect the gapped states in free 



In this paper, we studied the influence of different kinds 
of gaps in bilayer graphene in the two- and four-band 
models on longitudinal and transverse optical conductiv- 
ities paying special attention to gaps that break the time- 
reversal symmetry. The two-band model is valid for ener- 
gies E < 100 meV and the four-band model is applicable 
up to energies when continuum approximation is valid, 
i.e., for wave vectors ka -C 1 where a is the lattice con- 
stant. The upper bound 100 meV in the two-band model 
is due to neglecting the high-energy bands. Correspond- 
ing restrictions on gaps and on transverse voltages are 
related to limitations of these models. 

Starting from a low-energy two-band Hamiltonian, we 
derived a simple analytical expression (|14[) for the com- 
plex magneto-optical conductivity for two opposite circu- 
lar polarizations of light. It is verified that for zero values 
of the gaps, the strengths of the absorption lines satisfy 
the optical spectral weight conservation. When there is a 
layer asymmetry gap, the corresponding absorption peak 
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FIG. 7: (Color online) The calculated Faraday and the Kerr 
rotation angles for different experimental geometries shown 
in the inset. Optical conductivity of graphene was calculated 
using At = 1 meV and A = U = Ut = 0. 



Meanwhile the real part of the optical conductivity for 
the QAH state has a unique zero (see Fig. 4), the real 
part of the optical conductivity for the state with a su- 
perposition of LAF, QSH, and QVH gaps may have two 
zeros as a function of the energy of photon. Therefore, 
the observation of the optical Hall conductivity in zero 
magnetic field is a very effective probe of the ground state 
of bilaycr graphene. 

The time-reversal symmetry breaking states are ex- 
pected to be observed experimentally via optical polar- 
ization rotation either in the Faraday or Kerr effects. 
We analyzed a possibility of such experiments for a free 
standing graphene, graphene on a thick substrate and 
graphene on a double-layer substrate. In the last case 
the Faraday angle and especially the Kerr angle are no 
longer determined by the real part of the optical Hall 
conductivity only. Moreover, the sign of the Kerr an- 
gle is even inverted with respect to the case of a simple 
substrate. 

In the this paper, we calculated optical conductivities 
in a very simple approximation where quasiparticlc gaps 
were introduced into the Green's functions phenomeno- 
logically while the vertex corrections due to Coulomb 
and other interactions were neglected. Below the band 
gap the spectrum of the system may contain Coulomb 
bound electron-hole pairs which would reveal themselves 
as poles in the vertex function, hence as resonances in 
the absorption spectra. The spectrum of resonances de- 
pends on a type of dynamically generated gap and ex- 
perimental observation of these resonances could serve 
as another fingerprint for a given gapped ground state 
of bilayer graphene. Detailed study of such modes is be- 
yond of the scope of this paper and is postponed for the 
future. 



splits into two, while for the time-reversal symmetry gap, 
the position of the peak shifts, but it remains unsplit. 

The limit of zero magnetic field was analyzed for an 
arbitrary carrier density in the two-band approximation. 
We find that the necessary (but not sufficient) condi- 
tion for the optical Hall and dc Hall conductivities to 
remain finite is the presence of nonzero time-reversal 
symmetry breaking gap. In addition to the canonical 
time-reversal symmetry breaking QAH state considered 
in Ref. HU which provides Re a xy /(e 2 /h) = ±4 reflecting 
the presence of four topologically protected edge states 
for the QAH ground state of bilaycr graphene, we find 
another more sophisticated possibility of nonzero dc Hall 
conductivity Rea xy /(e 2 /h) = ±2 if time-reversal sym- 
metry breaking LAF gap, QSH gap, and QVH gap are 
present and satisfy a certain inequality. 

Using the full four-band model we derived analytic ex- 
pressions ([30j) and (|34|) for the optical Hall and longi- 
tudinal conductivities in a neutral bilayer graphene tak- 
ing into account the presence of four different gaps. We 
found that the optical Hall conductivity as a function of 
the energy of photon is strongly sensitive to the pres- 
ence of different time-reversal symmetry breaking states. 
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Appendix A: Diagonal conductivity spectral weight 

The analysis of the spectral weight proved to be useful 
both for theoretical (see e.g. Refs. I33T1391160I ) and exper- 
imental, Refl55l. studies of graphene. Here we consider 
the optical spectral weight that falls between O = and 
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O = Q m with n m a variable upper limit in the integral 

W(n m ) = / dnRea xx (n) (Al) 
Jo 

to verify that this weight is conserved irrespectively the 
value of the magnetic field. Setting A^ = in Eq. (fT6|) 
we arrive at the expression 



2c 



2 r 



| M |5(0) + ^(Q-2H) 



(A2) 



which could also be derived directly from the 2x2 
Hamiltonian^i for B = 0. Then for fi m > 2|/z| 



(A3) 



which is twice as much as the spectral weight for 
monolayer.— Obviously, Eqs. (|A2j) and (|A3j) are appli- 
cable only in the domain of validity of the 2x2 Hamil- 
tonian ([3]), e.g. for f2 < 71, but this is sufficient for our 
purposes. 

In the zero gap case, = Eq. (|14[) acquires a simple 
form 



To be specific, we have considered explicitly in this sec- 
tion only the case M < fi < M 2 . For \i &]M N , M N+1 [, 
N > 2, we can show that the missing spectral weight in 
the lines n < N is provided by the single intraband line 
at M N+1 - M N . 

Using a representation, similar to Eq. (|A5J) . for the 
conductivity in the limit T,T — > 0, but written for the 
case fj, g]Mtv, -Mjv+if under investigation, we obtain the 
optical spectral weight lost in units of (e 2 /h)ui c : 



t 



M, 



k+l 



N+l 



+ M 



N 



M, 



JV+l 



M, 



N 



(A8) 

The first term on the left hand side is the spectral weight 
from all lines that disappeared from k = to N — 2. The 
second term is due to the reduction in intensity by factor 
1/2 of the line at k = N — 1. The quantity on the right- 
hand side is the optical weight of the intraband line which 
has picked up all of the lost intensity. Any violation of 
the individual spectral weight of the lines would violate 
this conservation of the optical spectral weight. 



2e 2 h 3 



fe=0 



E 



nm 2 ! 4 

n F (XM k+1 ] 



n F (X'M k+2 ) 



A,A'=± 



X'Mi 



fc+2 



AM; 



k+l 



(A4) 



n T A'M fc+2 ± \M k+1 + 2iT 



Assuming for simplicity that Mq < \i < M2 we get in 
the limit T — > and T — > that the diagonal optical 
conductivity is 



2e 2 

Rea xx (n) = —u 2 c 
n 



5{tt - M k+2 - M fc+ i) 



k=0 



M k+2 + M k+1 



Accordingly for A^ = 0, the spectral weight equals 



(A5) 



w(n m ) 



2e 2 
IT ' 

2e 2 



: J2(k + i) 



k=0 
N 

^ Q 'VkT2 + Vk , 



M k+2 + Mfc+i 



(A6) 



where the maximal N is estimated from the condition 
n m = Mjv+2 + M N+ i w 2u c N. The sum over k in 
can be evaluated analytically, 



N 



E sr = lV(N + l)(N + 2)) 



N 



k=0 



sfk 



Vk 



(A7) 



Then, we again arrive at the result (|A3[) . obtained for 
B = confirming the conservation of the spectral weight. 
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